Systems genetics and bioinformatics analyses using ESR1-correlated genes identify potential candidates underlying female bone development.

Estrogen receptor α (ESR1) is involved in E2 signaling and plays a major role in postmenopausal bone loss. However, the molecular network underlying ESR1 has not been explored. We used systems genetics and bioinformatics to identify important genes associated with Esr1 in postmenopausal bone loss. We identified ~2300 Esr1-coexpressed genes in female BXD bone femur, functional analysis of which revealed 'osteoblast signaling' as the most enriched pathway. PPI network led to the identification of 25 'female bone candidates'. The gene-regulatory analysis revealed RUNX2 as a key TF. ANKRD1 and RUNX2 were significantly different between osteoporosis patients and healthy controls. Sp7, Col1a1 and Pth1r correlated with multiple femur bone phenotypes in BXD mice. miR-3121-3p targeted Csf1, Ankrd1, Sp7 and Runx2. ß-estradiol treatment markedly increased the expression of these candidates in mouse osteoblast. Our study revealed that Esr1-correlated genes Ankrd1, Runx2, Csf1 and Sp7 may play important roles in female bone development.

Human Breast Tissue Microbiota Reveals Unique Microbial Signatures that Correlate with Prognostic Features in Adult Ethiopian Women with Breast Cancer.

Breast cancer (BC) is the leading cause of cancer mortality among women in Ethiopia. Overall, women of African ancestry have the highest death toll due to BC compared to other racial/ethnic groups. The cause of the disparity in mortality is unclear. Recently, studies conducted in the United States and other high-income countries highlighted the role of microbial dysbiosis in BC initiation, tumor growth, and treatment outcome. However, the extent to which inter-individual differences in the makeup of microbiota are associated with clinical and histopathological outcomes in Ethiopian women has not been studied. The goal of our study was to profile the microbiome in breast tumor and normal adjacent to tumor (NAT) tissues of the same donor and to identify associations between microbial composition and abundance and clinicopathological factors in Ethiopian women with BC. We identified 14 microbiota genera in breast tumor tissues that were distinct from NAT tissues, of which Sphingobium, Anaerococcus, Corynebacterium, Delftia, and Enhydrobacter were most significantly decreased in breast tumors compared to NAT tissues. Several microbial genera significantly differed by clinicopathological factors in Ethiopian women with BC. Specifically, the genus Burkholderia more strongly correlated with aggressive triple negative (TNBC) and basal-like breast tumors. The genera Alkanindiges, Anoxybacillus, Leifsonia, and Exiguobacterium most strongly correlated with HER2-E tumors. Luminal A and luminal B tumors also correlated with Anoxybacillus but not as strongly as HER2-E tumors. A relatively higher abundance of the genus Citrobacter most significantly correlated with advanced-stage breast tumors compared to early-stage tumors. This is the first study to report an association between breast microbial dysbiosis and clinicopathological factors in Ethiopian women.

Exploring Early-Stage Retinal Neurodegeneration in Murine Pigmentary Glaucoma: Insights From Gene Networks and miRNA Regulation Analyses.

Purpose: Glaucoma is a group of heterogeneous optic neuropathies characterized by the progressive degeneration of retinal ganglion cells. However, the underlying mechanisms have not been understood completely. We aimed to elucidate the genetic network associated with the development of pigmentary glaucoma with DBA/2J (D2) mouse model of glaucoma and corresponding genetic control D2-Gpnmb (D2G) mice carrying the wild type (WT) Gpnmb allele. Methods: Retinas isolated from 13 D2 and 12 D2G mice were subdivided into 2 age groups: pre-onset (1-6 months: samples were collected at approximately 1-2, 2-4, and 5-6 months) and post-onset (7-15 months: samples were collected at approximately 7-9, 10-12, and 13-15 months) glaucoma were compared. Differential gene expression (DEG) analysis and gene-set enrichment analyses were performed. To identify micro-RNAs (miRNAs) that target Gpnmb, miRNA expression levels were correlated with time point matched mRNA expression levels. A weighted gene co-expression network analysis (WGCNA) was performed using the reference BXD mouse population. Quantitative real-time PCR (qRT-PCR) was used to validate Gpnmb and miRNA expression levels. Results: A total of 314 and 86 DEGs were identified in the pre-onset and post-onset glaucoma groups, respectively. DEGs in the pre-onset glaucoma group were associated with the crystallin gene family, whereas those in the post-onset group were related to innate immune system response. Of 1329 miRNAs predicted to target Gpnmb, 3 miRNAs (miR-125a-3p, miR-3076-5p, and miR-214-5p) were selected. A total of 47 genes demonstrated overlapping with the identified DEGs between D2 and D2G, segregated into their time-relevant stages. Gpnmb was significantly downregulated, whereas 2 out of 3 miRNAs were significantly upregulated (P < 0.05) in D2 mice at both 3-and 10-month time points. Conclusions: These findings suggest distinct gene-sets involved in pre-and post-glaucoma in the D2 mouse. We identified three miRNAs regulating Gpnmb in the development of murine pigmentary glaucoma.

Expression Levels of the Tnni3k Gene in the Heart Are Highly Associated with Cardiac and Glucose Metabolism-Related Phenotypes and Functional Pathways.

BACKGROUND: Troponin-I interacting kinase encoded by the TNNI3K gene is expressed in nuclei and Z-discs of cardiomyocytes. Mutations in TNNI3K were identified in patients with cardiac conduction diseases, arrhythmias, and cardiomyopathy. METHODS: We performed cardiac gene expression, whole genome sequencing (WGS), and cardiac function analysis in 40 strains of BXD recombinant inbred mice derived from C57BL/6J (B6) and DBA/2J (D2) strains. Expression quantitative trait loci (eQTLs) mapping and gene enrichment analysis was performed, followed by validation of candidate Tnni3k-regulatory genes. RESULTS: WGS identified compound splicing and missense T659I Tnni3k variants in the D2 parent and some BXD strains (D allele) and these strains had significantly lower Tnni3k expression than those carrying wild-type Tnni3k (B allele). Expression levels of Tnni3k significantly correlated with multiple cardiac (heart rate, wall thickness, PR duration, and T amplitude) and metabolic (glucose levels and insulin resistance) phenotypes in BXDs. A significant cis-eQTL on chromosome 3 was identified for the regulation of Tnni3k expression. Furthermore, Tnni3k-correlated genes were primarily involved in cardiac and glucose metabolism-related functions and pathways. Genes Nodal, Gnas, Nfkb1, Bmpr2, Bmp7, Smad7, Acvr1b, Acvr2b, Chrd, Tgfb3, Irs1, and Ppp1cb were differentially expressed between the B and D alleles. CONCLUSIONS: Compound splicing and T659I Tnni3k variants reduce cardiac Tnni3k expression and Tnni3k levels are associated with cardiac and glucose metabolism-related phenotypes.

Peer mentoring to support career advancement among underrepresented minority faculty in the programs to increase diversity among individuals engaged in health-related research (PRIDE).

Although mentoring is critical for career advancement, underrepresented minority (URM) faculty often lack access to mentoring opportunities. We sought to evaluate the impact of peer mentoring on career development success of URM early career faculty in the National Heart Lung and Blood Institute-sponsored, Programs to Increase Diversity Among Individuals Engaged in Health-Related Research-Functional and Translational Genomics of Blood Disorders (PRIDE-FTG). The outcome of peer mentoring was evaluated using the Mentoring Competency Assessment (MCA), a brief open-ended qualitative survey, and a semi-structured exit interview. Surveys were completed at baseline (Time 1), 6 months, and at the end of PRIDE-FTG participation (Time 2). The following results were obtained. Between Time 1 and Time 2, mentees' self-assessment scores increased for the MCA (p < 0.01) with significant increases in effective communication (p < 0.001), aligning expectations (p < 0.05), assessing understanding (p < 0.01), and addressing diversity (p < 0.002). Mentees rated their peer mentors higher in the MCA with significant differences noted for promoting development (p < 0.027). These data suggest that PRIDE-FTG peer mentoring approaches successfully improved MCA competencies among URM junior faculty participants with faculty ranking peer mentors higher than themselves. Among URM faculty, peer mentoring initiatives should be investigated as a key strategy to support early career scholar development.

Weighted gene co-expression network analysis identifies key hub genes and pathways in acute myeloid leukemia.

Introduction: Acute myeloid leukemia (AML) is the most common type of leukemia in adults. However, there is a gap in understanding the molecular basis of the disease, partly because key genes associated with AML have not been extensively explored. In the current study, we aimed to identify genes that have strong association with AML based on a cross-species integrative approach. Methods: We used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify co-expressed gene modules significantly correlated with human AML, and further selected the genes exhibiting a significant difference in expression between AML and healthy mouse. Protein-protein interactions, transcription factors, gene function, genetic regulation, and coding sequence variants were integrated to identify key hub genes in AML. Results: The cross-species approach identified a total of 412 genes associated with both human and mouse AML. Enrichment analysis confirmed an association of these genes with hematopoietic and immune-related functions, phenotypes, processes, and pathways. Further, the integrated analysis approach identified a set of important module genes including Nfe2, Trim27, Mef2c, Ets1, Tal1, Foxo1, and Gata1 in AML. Six of these genes (except ETS1) showed significant differential expression between human AML and healthy samples in an independent microarray dataset. All of these genes are known to be involved in immune/hematopoietic functions, and in transcriptional regulation. In addition, Nfe2, Trim27, Mef2c, and Ets1 harbor coding sequence variants, whereas Nfe2 and Trim27 are cis-regulated, making them attractive candidates for validation. Furthermore, subtype-specific analysis of the hub genes in human AML indicated high expression of NFE2 across all the subtypes (M0 through M7) and enriched expression of ETS1, LEF1, GATA1, and TAL1 in M6 and M7 subtypes. A significant correlation between methylation status and expression level was observed for most of these genes in AML patients. Conclusion: Findings from the current study highlight the importance of our cross-species approach in the identification of multiple key candidate genes in AML, which can be further studied to explore their detailed role in leukemia/AML.

Cardiac copper content and its relationship with heart physiology: Insights based on quantitative genetic and functional analyses using BXD family mice.

Background: Copper (Cu) is essential for the functioning of various enzymes involved in important cellular and physiological processes. Although critical for normal cardiac function, excessive accumulation, or deficiency of Cu in the myocardium is detrimental to the heart. Fluctuations in cardiac Cu content have been shown to cause cardiac pathologies and imbalance in systemic Cu metabolism. However, the genetic basis underlying cardiac Cu levels and their effects on heart traits remain to be understood. Representing the largest murine genetic reference population, BXD strains have been widely used to explore genotype-phenotype associations and identify quantitative trait loci (QTL) and candidate genes. Methods: Cardiac Cu concentration and heart function in BXD strains were measured, followed by QTL mapping. The candidate genes modulating Cu homeostasis in mice hearts were identified using a multi-criteria scoring/filtering approach. Results: Significant correlations were identified between cardiac Cu concentration and left ventricular (LV) internal diameter and volumes at end-diastole and end-systole, demonstrating that the BXDs with higher cardiac Cu levels have larger LV chamber. Conversely, cardiac Cu levels negatively correlated with LV posterior wall thickness, suggesting that lower Cu concentration in the heart is associated with LV hypertrophy. Genetic mapping identified six QTLs containing a total of 217 genes, which were further narrowed down to 21 genes that showed a significant association with cardiac Cu content in mice. Among those, Prex1 and Irx3 are the strongest candidates involved in cardiac Cu modulation. Conclusion: Cardiac Cu level is significantly correlated with heart chamber size and hypertrophy phenotypes in BXD mice, while being regulated by multiple genes in several QTLs. Prex1 and Irx3 may be involved in modulating Cu metabolism and its downstream effects and warrant further experimental and functional validations.

MicroRNA29B induces fetal hemoglobin via inhibition of the HBG repressor protein MYB in vitro and in humanized sickle cell mice.

Introduction: Therapeutic strategies aimed at reactivating HBG gene transcription and fetal hemoglobin (HbF) synthesis remain the most effective strategy to ameliorate the clinical symptoms of sickle cell disease (SCD). We previously identified microRNA29B (MIR29B) as a novel HbF inducer via targeting enzymes involved in DNA methylation. We provided further evidence that the introduction of MIR29B into KU812 leukemia cells significantly reduced MYB protein expression. Therefore, the aim of this study was to determine the extent to which MIR29B mediates HbF induction via targeting MYB in KU812 leukemia cells and human primary erythroid progenitors and to investigate the role of MIR29B in HbF induction in vivo in the humanized Townes SCD mouse model. Materials and methods: Human KU812 were cultured and normal CD34 cells (n = 3) were differentiated using a two-phase erythropoiesis culturing system and transfected with MIR29B (50 and 100 nM) mimic or Scrambled (Scr) control in vitro. A luciferase reporter plasmid overexpressing MYB was transfected into KU812 cells. Luciferase activity was quantified after 48 h. Gene expression was determined by quantitative real-time PCR. In vivo studies were conducted using Townes SCD mice (6 per group) treated with MIR29B (2, 3, and 4 mg/kg/day) or Scr control by 28-day continuous infusion using subcutaneous mini osmotic pumps. Blood samples were collected and processed for complete blood count (CBC) with differential and reticulocytes at weeks 0, 2, and 4. Flow cytometry was used to measure the percentage of HbF-positive cells. Results: In silico analysis predicted complementary base-pairing between MIR29B and the 3'-untranslated region (UTR) of MYB. Overexpression of MIR29B significantly reduced MYB mRNA and protein expression in KU812 cells and erythroid progenitors. Using a luciferase reporter vector that contained the full-length MYB 3'-UTR, we observed a significant reduction in luciferase activity among KU812 cells that co-expressed MIR29B and the full-length MYB 3'-UTR as compared to cells that only expressed MYB 3'-UTR. We confirmed the inhibitory effect of a plasmid engineered to overexpress MYB on HBG activation and HbF induction in both KU812 cells and human primary erythroid progenitors. Co-expression of MIR29B and MYB in both cell types further demonstrated the inhibitory effect of MIR29B on MYB expression, resulting in HBG reactivation by real-time PCR, Western blot, and flow cytometry analysis. Finally, we confirmed the ability of MIR29B to reduce sickling and induce HbF by decreasing expression of MYB and DNMT3 gene expression in the humanized Townes sickle cell mouse model. Discussion: Our findings support the ability of MIR29B to induce HbF in vivo in Townes sickle cell mice. This is the first study to provide evidence of the ability of MIR29B to modulate HBG transcription by MYB gene silencing in vivo. Our research highlights a novel MIR-based epigenetic approach to induce HbF supporting the discovery of new drugs to expand treatment options for SCD.

Characterization of the Metabolome of Breast Tissues from Non-Hispanic Black and Non-Hispanic White Women Reveals Correlations between Microbial Dysbiosis and Enhanced Lipid Metabolism Pathways in Triple-Negative Breast Tumors.

Triple-negative breast cancer (TNBC) is an aggressive form of breast cancer that is non-responsive to hormonal therapies and disproportionately impact women of African ancestry. We previously showed that TN breast tumors have a distinct microbial signature that differs from less aggressive breast tumor subtypes and normal breast tissues. However, it is unknown whether these differences in breast tumor microbiota may be driven by alterations in microbial metabolites, leading to potentially protective or pathogenic consequences. The goal of this global metabolomic profiling study was to investigate alterations in microbial metabolism pathways in normal and breast tumor tissues, including TNBC, of non-Hispanic black (NHB) and non-Hispanic white (NHW) women. In this study, we profiled the microbiome (16S rRNA) from breast tumor tissues and analyzed 984 metabolites from a total of 51 NHB and NHW women. Breast tumor tissues were collected from 15 patients with TNBC, 12 patients with less aggressive luminal A-type (Luminal) breast cancer, and 24 healthy controls for comparison using UHPLC-tandem mass spectrometry. Principal component analysis and hierarchical clustering of the global metabolomic profiling data revealed separation between metabolic signatures of normal and breast tumor tissues. Random forest analysis revealed a unique biochemical signature associated with elevated lipid metabolites and lower levels of microbial-derived metabolites important in controlling inflammation and immune responses in breast tumor tissues. Significant relationships between the breast microbiome and the metabolome, particularly lipid metabolism, were observed in TNBC tissues. Further investigations to determine whether alterations in sphingolipid, phospholipid, ceramide, amino acid, and energy metabolism pathways modulate Fusobacterium and Tenericutes abundance and composition to alter host metabolism in TNBC are necessary to help us understand the risk and underlying mechanisms and to identify potential microbial-based targets.

Targeting Genetic Modifiers of HBG Gene Expression in Sickle Cell Disease: The miRNA Option.

Sickle cell disease (SCD) is one of the most common inherited hemoglobinopathy disorders that affects millions of people worldwide. Reactivation of HBG (HBG1, HBG2) gene expression and induction of fetal hemoglobin (HbF) is an important therapeutic strategy for ameliorating the clinical symptoms and severity of SCD. Hydroxyurea is the only US FDA-approved drug with proven efficacy to induce HbF in SCD patients, yet serious complications have been associated with its use. Over the last three decades, numerous additional pharmacological agents that reactivate HBG transcription in vitro have been investigated, but few have proceeded to FDA approval, with the exception of arginine butyrate and decitabine; however, neither drug met the requirements for routine clinical use due to difficulties with oral delivery and inability to achieve therapeutic levels. Thus, novel approaches that produce sufficient efficacy, specificity, and sustainable HbF induction with low adverse effects are desirable. More recently, microRNAs (miRNAs) have gained attention for their diagnostic and therapeutic potential to treat various diseases ranging from cancer to Alzheimer's disease via targeting oncogenes and their gene products. Thus, it is plausible that miRNAs that target HBG regulatory genes may be useful for inducing HbF as a treatment for SCD. Our laboratory and others have documented the association of miRNAs with HBG activation or suppression via silencing transcriptional repressors and activators, respectively, of HBG expression. Herein, we review progress made in understanding molecular mechanisms of miRNA-mediated HBG regulation and discuss the extent to which molecular targets of HBG might be suitable prospects for development of SCD clinical therapy. Lastly, we discuss challenges with the application of miRNA delivery in vivo and provide potential strategies for overcoming barriers in the future.

A systems genetics approach delineates the role of Bcl2 in leukemia pathogenesis.

Leukemia is a group of malignancies of the blood forming tissues, and is characterized by the uncontrolled proliferation of blood cells. In the United States, it accounts for approximately 3.5% and 4% of all cancer-related incidences and mortalities, respectively. The current study aimed to explore the role of Bcl2 and associated genes in leukemia pathogenesis using a systems genetics approach. The transcriptome data from BXD Recombinant Inbred (RI) mice was analyzed to identify the expression of Bcl2 in myeloid cells. eQTL mapping was performed to select the potential chromosomal region and subsequently identify the candidate gene modulating the expression of Bcl2. Furthermore, gene enrichment and protein-protein interaction (PPI) analyses of the Bcl2-coexpressed genes were performed to demonstrate the role of Bcl2 in leukemia pathogenesis. The Bcl2-coexpressed genes were found to be enriched in various hematopoietic system related functions, and multiple pathways related to signaling, immune response, and cancer. The PPI network analysis demonstrated direct interaction of hematopoietic function related genes, such as Bag3, Bak1, Bcl2l11, Bmf, Mapk9, Myc, Ppp2r5c, and Ppp3ca with Bcl2. The eQTL mapping identified a 4.5 Mb genomic region on chromosome 11, potentially regulating the expression of Bcl2. A multi-criteria filtering process identified Top2a, among the genes located in the mapped locus, as the best candidate upstream regulator for Bcl2 expression variation. Hence, the current study provides better insights into the role of Bcl2 in leukemia pathogenesis and demonstrates the significance of our approach in gaining new knowledge on leukemia. Furthermore, our findings from the PPI network analysis and eQTL mapping provide supporting evidence of leukemia-associated genes, which can be further explored for their functional importance in leukemia. DATA AVAILABILITY: The myeloid cell transcriptomic data of the BXD mice used in this study can be accessed through our GeneNetwork (http://www.genenetwork.org) with the accession number of GN144.

Exploring epigenetic and microRNA approaches for γ-globin gene regulation.

Therapeutic interventions aimed at inducing fetal hemoglobin and reducing the concentration of sickle hemoglobin is an effective approach to ameliorating acute and chronic complications of sickle cell disease, exemplified by the long-term use of hydroxyurea. However, there remains an unmet need for the development of additional safe and effective drugs for single agent or combination therapy for individuals with ß-hemoglobinopathies. Regulation of the γ-globin to ß-globin switch is achieved by chromatin remodeling at the HBB locus on chromosome 11 and interactions of major DNA binding proteins, such as KLF1 and BCL11A in the proximal promoters of the globin genes. Experimental evidence also supports a role of epigenetic modifications including DNA methylation, histone acetylation/methylation, and microRNA expression in γ-globin gene silencing during development. In this review, we will critically evaluate the role of epigenetic mechanisms in γ-globin gene regulation and discuss data generated in tissue culture, pre-clinical animal models, and clinical trials to support drug development to date. The question remains whether modulation of epigenetic pathways will produce sufficient efficacy and specificity for fetal hemoglobin induction and to what extent targeting these pathways form the basis of prospects for clinical therapy.

Sickle cell disease: progress towards combination drug therapy.

Dr. John Herrick described the first clinical case of sickle cell anaemia (SCA) in the United States in 1910. Subsequently, four decades later, Ingram and colleagues characterized the A to T substitution in DNA producing the GAG to GTG codon and replacement of glutamic acid with valine in the sixth position of the ßS -globin chain. The establishment of Comprehensive Sickle Cell Centers in the United States in the 1970s was an important milestone in the development of treatment strategies and describing the natural history of sickle cell disease (SCD) comprised of genotypes including homozygous haemoglobin SS (HbSS), HbSß0 thalassaemia, HbSC and HbSß+ thalassaemia, among others. Early drug studies demonstrating effective treatments of HbSS and HbSß0 thalassaemia, stimulated clinical trials to develop disease-specific therapies to induce fetal haemoglobin due to its ability to block HbS polymerization. Subsequently, hydroxycarbamide proved efficacious in adults with SCA and was Food and Drug Administration (FDA)-approved in 1998. After two decades of hydroxycarbamide use for SCD, there continues to be limited clinical acceptance of this chemotherapy drug, providing the impetus for investigators and pharmaceutical companies to develop non-chemotherapy agents. Investigative efforts to determine the role of events downstream of deoxy-HbS polymerization, such as endothelial cell activation, cellular adhesion, chronic inflammation, intravascular haemolysis and nitric oxide scavenging, have expanded drug targets which reverse the pathophysiology of SCD. After two decades of slow progress in the field, since 2018 three new drugs were FDA-approved for SCA, but research efforts to develop treatments continue. Currently over 30 treatment intervention trials are in progress to investigate a wide range of agents acting by complementary mechanisms, providing the rationale for ushering in the age of effective and safe combination drug therapy for SCD. Parallel efforts to develop curative therapies using haematopoietic stem cell transplant and gene therapy provide individuals with SCD multiple treatment options. We will discuss progress made towards drug development and potential combination drug therapy for SCD with the standard of care hydroxycarbamide.

The Genetic Dissection of Ace2 Expression Variation in the Heart of Murine Genetic Reference Population.

Background: A high inflammatory and cytokine burden that induces vascular inflammation, myocarditis, cardiac arrhythmias, and myocardial injury is associated with a lethal outcome in COVID-19. The SARS-CoV-2 virus utilizes the ACE2 receptor for cell entry in a similar way to SARS-CoV. This study investigates the regulation, gene network, and associated pathways of ACE2 that may be involved in inflammatory and cardiovascular complications of COVID-19. Methods: Cardiovascular traits were determined in the one of the largest mouse genetic reference populations: BXD recombinant inbred strains using blood pressure, electrocardiography, and echocardiography measurements. Expression quantitative trait locus (eQTL) mapping, genetic correlation, and functional enrichment analysis were used to identify Ace2 regulation, gene pathway, and co-expression networks. Results: A wide range of variation was found in expression of Ace2 among the BXD strains. Levels of Ace2 expression are negatively correlated with cardiovascular traits, including systolic and diastolic blood pressure and P wave duration and amplitude. Ace2 co-expressed genes are significantly involved in cardiac- and inflammatory-related pathways. The eQTL mapping revealed that Cyld is a candidate upstream regulator for Ace2. Moreover, the protein-protein interaction (PPI) network analysis inferred several potential key regulators (Cul3, Atf2, Vcp, Jun, Ppp1cc, Npm1, Mapk8, Set, Dlg1, Mapk14, and Hspa1b) for Ace2 co-expressed genes in the heart. Conclusions: Ace2 is associated with blood pressure, atrial morphology, and sinoatrial conduction in BXD mice. Ace2 co-varies with Atf2, Cyld, Jun, Mapk8, and Mapk14 and is enriched in the RAS, TGFß, TNFα, and p38α signaling pathways, involved in inflammation and cardiac damage. We suggest that all these novel Ace2-associated genes and pathways may be targeted for preventive, diagnostic, and therapeutic purposes in cardiovascular damage in patients with systemic inflammation, including COVID-19 patients.

Exploring the Role of Chemokine Receptor 6 (Ccr6) in the BXD Mouse Model of Gulf War Illness.

Gulf War illness (GWI) is a chronic and multi-symptomatic disorder with persistent neuroimmune symptomatology. Chemokine receptor 6 (CCR6) has been shown to be involved in several inflammation disorders in humans. However, the causative relationship between CCR6 and neuroinflammation in GWI has not yet been investigated. By using RNA-seq data of prefrontal cortex (PFC) from 31 C57BL/6J X DBA/2J (BXD) recombinant inbred (RI) mouse strains and their parental strains under three chemical treatment groups - saline control (CTL), diisopropylfluorophosphate (DFP), and corticosterone combined with diisopropylfluorophosphate (CORT+DFP), we identified Ccr6 as a candidate gene underlying individual differences in susceptibility to GWI. The Ccr6 gene is cis-regulated and its expression is significantly correlated with CORT+DFP treatment. Its mean transcript abundance in PFC of BXD mice decreased 1.6-fold (p < 0.0001) in the CORT+DFP group. The response of Ccr6 to CORT+DFP is also significantly different (p < 0.0001) between the parental strains, suggesting Ccr6 is affected by both host genetic background and chemical treatments. Pearson product-moment correlation analysis revealed 1473 Ccr6-correlated genes (p < 0.05). Enrichment of these genes was seen in the immune, inflammation, cytokine, and neurological related categories. In addition, we also found five central nervous system-related phenotypes and fecal corticosterone concentration have significant correlation (p < 0.05) with expression of Ccr6 in the PFC. We further established a protein-protein interaction subnetwork for the Ccr6-correlated genes, which provides an insight on the interaction of G protein-coupled receptors, kallikrein-kinin system and neuroactive ligand-receptors. This analysis likely defines the heterogeneity and complexity of GWI. Therefore, our results suggest that Ccr6 is one of promising GWI biomarkers.

An Age-Related Hearing Protection Locus on Chromosome 16 of BXD Strain Mice.

Inbred mouse models are widely used to study age-related hearing loss (AHL). Many genes associated with AHL have been mapped in a variety of strains. However, little is known about gene variants that have the converse function-protective genes that confer strong resistance to hearing loss. Previously, we reported that C57BL/6J (B6) and DBA/2J (D2) strains share a common hearing loss allele in Cdh23. The cadherin 23 (Cdh23) gene is a key contributor to early-onset hearing loss in humans. In this study, we tested hearing across a large family of 54 BXD strains generated from B6 to D2 crosses. Five of 54 strains maintain the normal threshold (20 dB SPL) even at 2 years old-an age at which both parental strains are essentially deaf. Further analyses revealed an age-related hearing protection (ahp) locus on chromosome 16 (Chr 16) at 57~76 Mb with a maximum LOD of 5.7. A small number of BXD strains at 2 years with good hearing correspond roughly to the percentage of humans who have good hearing at 90 years old. Further studies to define candidate genes in the ahp locus and related molecular mechanisms involved in age-related resilience or resistance to AHL are warranted.

Genome-wide transcriptome architecture in a mouse model of Gulf War Illness.

Gulf War Illness (GWI) is thought to be a chronic neuroimmune disorder caused by in-theater exposure during the 1990-1991 Gulf War. There is a consensus that the illness is caused by exposure to insecticides and nerve agent toxicants. However, the heterogeneity in both development of disease and clinical outcomes strongly suggests a genetic contribution. Here, we modeled GWI in 30 BXD recombinant inbred mouse strains with a combined treatment of corticosterone (CORT) and diisopropyl fluorophosphate (DFP). We quantified transcriptomes from 409 prefrontal cortex samples. Compared to the untreated and DFP treated controls, the combined treatment significantly activated pathways such as cytokine-cytokine receptor interaction and TNF signaling pathway. Protein-protein interaction analysis defined 6 subnetworks for CORT + DFP, with the key regulators being Cxcl1, Il6, Ccnb1, Tnf, Agt, and Itgam. We also identified 21 differentially expressed genes having significant QTLs related to CORT + DFP, but without evidence for untreated and DFP treated controls, suggesting regions of the genome specifically involved in the response to CORT + DFP. We identified Adamts9 as a potential contributor to response to CORT + DFP and found links to symptoms of GWI. Furthermore, we observed a significant effect of CORT + DFP treatment on the relative proportion of myelinating oligodendrocytes, with a QTL on Chromosome 5. We highlight three candidates, Magi2, Sema3c, and Gnai1, based on their high expression in the brain and oligodendrocyte. In summary, our results show significant genetic effects of the CORT + DFP treatment, which mirrors gene and protein expression changes seen in GWI sufferers, providing insight into the disease and a testbed for future interventions.

Pex3 is involved in the genetic regulation of Nr3c2 expression in the amygdala of mice.

The mineralocorticoid receptor (Nr3c2) has received increased attention as an important stress-related gene. Here, we sought to uncover candidate genes regulating the expression of Nr3c2. Using a genetical genomics approach, we identified a significant trans-regulated expression quantitative trait locus (eQTL) at Chromosome 10 for Nr3c2 expression in the amygdala of BXD RI strains. We then examined genes upstream of the eQTL to identify likely regulatory candidates of Nr3c2 expression. Pex3 (peroxisomal) expression was highly correlated with that of Nr3c2, had a significant cis-regulated eQTL that mapped to the Nr3c2 eQTL region and thus emerged as the most likely regulatory candidate of Nr3c2 expression. In vitro studies showed that silencing of Pex3 by siRNA decreased Nr3c2 expression in HEK293T and SHSY5 cell lines while overexpression increased Nr3c2 expression. A relationship between the expression of these two genes was further supported by our observations that expression levels of Pex3 and Nr3c2 decreased in the amygdala of mice exposed to chronic unpredictable stress. Our findings provide insight into the genetic regulation of Nr3c2 expression and suggest a new role for Pex3 in stress responses. Future characterization of Pex3's role in the regulation of Nr3c2 expression and the pathways involved may lead to a better understanding of stress responses and risk for stress-related pathology.

Prevalence of the UGT1A1*28 promoter polymorphism and breast cancer risk among African American women in Memphis, TN.

Inherited variations in UDP-glucuronosyltransferase 1A1 (UGT1A1) are associated with an increased breast cancer risk in women of African ancestry. The UGT1A1*28 promoter polymorphism is characterized by the presence of 7 TA repeats in the TATA box sequence and results in reduced UGT1A1 gene expression and enzymatic activity. In this study, we investigated associations between the UGT1A1*28 polymorphism and breast cancer risk among African American (AA) women in Memphis, Tennessee, a city with increased breast cancer mortality rates among AA women. Saliva was collected from 352 AA women, including breast cancer cases (n=82) and controls (n=270) between June 2016 to June 2017. DNA was isolated and sequenced for the UGT1A1*28 polymorphism. The odds ratio for cases with the low UGT1A1 activity alleles (TA)7/8 repeat genotypes versus 5/5, 5/6, and 6/6 genotypes was 1.46 [95% CI, 0.65-3.31; P = 0.36] in premenopausal women and 1.10 (95% CI, 0.52-2.38; P = 0.79) in postmenopausal women. Further analysis of TCGA RNA-seq data showed that UGT1A1 mRNA was significantly lower among estrogen receptor (ER)-negative breast cancers from AA as compared to non-Hispanic white women with ER-negative breast cancer. Larger epidemiological studies are needed to determine the functional consequence of the UGT1A1*28 polymorphism on breast cancer risk in AA women.

Distinct microbial communities that differ by race, stage, or breast-tumor subtype in breast tissues of non-Hispanic Black and non-Hispanic White women.

Growing evidence highlights an association between an imbalance in the composition and abundance of bacteria in the breast tissue (referred as microbial dysbiosis) and breast cancer in women. However, studies on the breast tissue microbiome have not been conducted in non-Hispanic Black (NHB) women. We investigated normal and breast cancer tissue microbiota from NHB and non-Hispanic White (NHW) women to identify distinct microbial signatures by race, stage, or tumor subtype. Using 16S rRNA gene sequencing, we observed that phylum Proteobacteria was most abundant in normal (n = 8), normal adjacent to tumor (normal pairs, n = 11), and breast tumors from NHB and NHW women (n = 64), with fewer Firmicutes, Bacteroidetes, and Actinobacteria. Breast tissues from NHB women had a higher abundance of genus Ralstonia compared to NHW tumors, which could explain a portion of the breast cancer racial disparities. Analysis of tumor subtype revealed enrichment of family Streptococcaceae in TNBC. A higher abundance of genus Bosea (phylum Proteobacteria) increased with stage. This is the first study to identify racial differences in the breast tissue microbiota between NHB and NHW women. Further studies on the breast cancer microbiome are necessary to help us understand risk, underlying mechanisms, and identify potential microbial targets.